Dynamical mean field theory of small polaron transport 
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We present a unified view of the transport properties of small-polarons in the Holstein model at 
low carrier densities, based on the Dynamical Mean Field Theory. The nonperturbative nature of 
the approach allows us to study the crossover from classical activated motion at high temperatures to 
coherent motion at low temperatures. Large quantitative discrepancies from the standard polaronic 
formulae are found. The scaling properties of the resistivity are analysed, and a simple interpolation 
formula is proposed in the nonadiabatic regime. 
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In the common wisdom, polaronic transport in solids 
is synonymous of activated conductivity. When elec- 
trons are strongly coupled to bosonic degrees of freedom 
(phonons, excitons, etc.), self-trapped states are formed. 
If the size of the polaron is comparable with the lattice 
spacing, the motion is then dominated by hopping pro- 
cesses: the particle has to overcome a potential barrier 
A and loses its quantum coherence at each hop, giving 
rise to an activated law of the form: 



P = Poe 



A/k B T 



(1) 



where the prefactor po is weakly temperature dependent. 
Contrary to what happens in semiconductors, the acti- 
vated behaviour here is not related to the number of ther- 
mally excited carriers, but rather to the mobility of the 
individual carriers. 
The Holstein modeli^ 



H = -f y^(c|cj + c\c,j) + huj a\ai 



was introduced in the late fifties to study such a be- 
haviour, as was being measured in some transition metal 
oxides. In this model, tight-binding electrons interact 
locally with molecular deformations, whose natural vi- 
bration frequency is luq [t is the hopping parameter, g 
the electron-phonon coupling constant). Although this 
is a rather crude idealization of a real solid, the Hol- 
stein model captures the essential physical phenomena 
involved in small-polaron transport. The situation is 
in fact more complex than indicated by the simple for- 
mula Q), and is summarized in several reviews^. Three 
regimes of temperature can be identified. 

At low temperatures, the polarons behave as heavy 
particles in a band of renormalized width W, and 
are weakly scattered by phonons. For small enough 
W, all the states in the band are equally popu- 
lated, leading to a "metallic-like" resistivity^^ p oc 
(T/W) exp(— TiuJo/kBT). The exponential law comes 
from the thermal occupation of the optical phonons, 
which are assumed to be the main source of scattering. 



Upon increasing the temperature, the mean free path 
is rapidly reduced until the picture of coherent motion 
breaks down, typically around some fraction of hujo/kB- 
Hopping motion then becomes more favorable, leading 
to an activated behavior of the form JIJ. The crossover 
from coherent to hopping motion is thus characterised 
by a maximum in the resistivity. Eventually, at tem- 
peratures higher than the activation energy, the polaron 
states are thermally dissociated and the residual elec- 
trons are strongly scattered by thermal phonons. In this 
case, the equipartition principle leads tc- 6 p ~ T 3 / 2 . 

On the experimental side, the largest amount of work 
has been devoted to the activated regime, which is often 
observed around room temperature, and resistivities of 
the form Q have been measured in a variety of narrow- 
band solids. However, strong deviations from pure Arrhe- 
nius behaviour are often reported^, possibly indicating 
the onset of the coherent transport regime. In some cases, 
the low temperature exponential law described above has 
also been identified unambiguously^. 

The main purpose of this work is to shed some light 
on the crossover from activated to coherent transport, 
for which a reliable theoretical description is still lacking. 
To do this, we calculate the resistivity of Holstein small- 
polarons in the framework of the Dynamical Mean Field 
Theory (DMFT). This approximation is suitable when- 
ever the physics is ruled by local phenomena, as is the 
case in the present problem, where it allows to take into 
account the quantum nature of the phonons (luq ^ 0) and 
the finite bandwidth effects (t ^ 0) on the same footing. 
Since the theory does not require any "small parameter" , 
it is able to go beyond the traditional approaches usually 
applied to the problem and gives reliable results in the 
regime ksT ~ Too® of interest here. Moreover, it yields 
a unified view of the different regimes of polaronic trans- 
port, acting as a testing ground of the validity of previous 
approaches 

The DMFT solution of the Holstein model for a sin- 
gle polaron was presented in ref^S, where an analyti- 
cal expression for the spectral function A e (u) was given 
in terms of a continued fraction. The polaron forma- 
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tion at zero temperature can be described by introducing 
two independent dimensionless parameters. The first is 
the adiabaticity ratio 7 = ujq/D (D is the unrenormal- 
ized half bandwidth) according to which an adiabatic 
(7 -C f) and nonadiabatic regime (7 > 1) can be de- 
fined. The mechanism of polaron formation is funda- 
mentally different in the two regimes, leading to differ- 
ent definitions of the dimensionless electron-phonon cou- 
pling. Being Ep — g 2 /uiQ the energy of a polaron on a 
single lattice site, a well defined polaronic state is formed 
for large A = Ep/D in the adiabatic case, and for large 
a 2 = Ep/ujQ in the non adiabatic cas o 10 ! 11 . 

The corresponding transport properties can be calcu- 
lated through the appropriate Kubo formula, which re- 
lates them to the current-current correlation function of 
the system at equilibrium. In DMFT, due to the absence 
of vertex corrections^ the latter is fully determined by 
the spectral function A £ (v), which is known exactly in 
the limit of vanishing density (single polaron problem) 10 . 
The resistivity at low (but finite) density can then be de- 
rived through an expansion in the fugacity^, yielding 



P(T) 



k B T j deN e j dve- v l T A e (y) 



(2) 



x(n J deN t cj} e J dve~"l T [A t (v)Y 

In the above formula, the denominator is the current- 
current correlation function, and the numerator is pro- 
portional the carrier concentration x, which was ex- 
plicitely taken out (the exponential weights in the in- 
tegrals reflect the Boltzmann nature of the carriers) . N e 
and <j) t are respectively the density of states and the 
current vertex of the periodic latticed The constant 
£ = e 2 a 2 /hv has the dimensions of resistivity, a being 
the lattice spacing, v the volume of the unit cell. In 
the following, we shall implicitely report the results for 
the dimensionless product pxQ, which is inversely propor- 
tional to the drift mobility, and assume h — ks = 1. 
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FIG. 1: Resistivity vs. temperature at A = 1.5, for different 
values of the adiabaticity parameter. Short arrows mark the 
temperature T = ujq/2 below which phonon quantum effects 
acquire importance. The long upward arrow is the gap T = A. 

The results for the resistivity are illustrated in Fig. 
n for fixed A = 1.5 in the adiabatic case. The three 



regimes discussed in the introduction can be clearly iden- 
tified — the resistivity first rises exponentially (coherent 
regime) , then decreases exponentially (activated regime) , 
and eventually increases again as a power law (residual 
scattering regime). This is true for all the data sets ex- 
cept at 7 = 0.4, where the polaron formation has shifted 
to higher values of A, as expected when moving away 
from the adiabatic limit ^Lii 

Let us focus on the activated regime, wo ^ T < A, 
where the polaron transport is dominated by incoherent 
hopping processes. In the adiabatic limit 7 — > 0, the 
problem is generally studied within a simplified "two- 
site" molecular modeli^, where the (classical) lattice 
degrees of freedom are seen to move "adiabatically" in 
the double- well energy curve determined by the electronic 
(bonding) state. At each jump, the system has to over- 
come a barrier A = Ep/2 — t, leading to an Arrhcnius 
type behavior 3 



p = 2{T/u Q ) exp[A/T] 



(3) 



This semi-classical description holds provided that the 
transitions to higher (antibonding) electronic states can 
be neglected, which corresponds to^ 



D 2 



r\i 



LUoVZEpT 



= [2A 7 3 (I> )] 



-1/2 



> 1 



(4) 



Note that this does not coincide with the usual adiabatic- 
ity condition 7 <§; 1 relevant for polaron formation. 
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FIG. 2: Arrhenius plots of resistivity divided by T/ujq for 
fixed 7 = 0.2, at various A. The DMFT data are compared 
with the semi-classical formula © with A = [Ep(X) — D]/2 
(straight lines) — see text. 

To illustrate the accuracy of the semi-classical predic- 
tion, we show in Fig. [3 Arrhenius plots of the resistiv- 
ity at fixed 7 = 0.2, varying the coupling strength A. 
First of all, our results indicate that the correct gen- 
eralization of the "two-site" result (0 to infinite lat- 
tices is obtained by letting A = [Ep(X) — D]/2, where 
Ep(X) = DX + D/ (8A) + • • • is the adiabatic polaron en- 
ergy, calculated e.g. in ref>i°« [see Fig. 2, full lines — the 
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slight discrepancy at the highest values of A is related to 
the breakdown of the condition Q]. This suggests that 
the activated behavior arises from the thermal promotion 
from the ground state to the electron continuum, which 
differs from the Landau-Zener mechanism involved in the 
two-site model. In particular, the reduction of the acti- 
vation gap by finite bandwidth effects is much stronger 
in the present case. 

When the temperature is lowered below T w lu , the 
quantum nature of the phonons can no longer be ne- 
glected. The lattice zero point fluctuations are then ex- 
pected to induce derealization of the trapped carriers, 
resulting in an enhancement of the polaron mobility. This 
phenomenon is a precursor of the coherent regime, and 
has a sizeable influence on the transport properties in a 
wide range of temperatures between the resistivity max- 
imum Tb and the phonon frequency loq, which is in prin- 
ciple experimentally accessible. The enhancement of the 
mobility is signaled by a marked downturn from the Ar- 
rhenius behavior (see the right hand side of Fig. 2), and 
takes place in the whole polaronic regime A 1 (as we 
shall see below, this behavior is quite general and is not 
restricted to the adiabatic case). 
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FIG. 3: Resistivity vs. temperature at a 2 = 10. The data 
points correspond respectively to 7 =4 (+), 2 (squares), 1 (x), 
0.5 (open circles), 0.2 (triangles), 0.15 (filled circles), 0.1 (open 
triangles) and were divided by 7 to evidence the nonadiabatic 
scaling property @. Inset: Arrhenius plot of p at 7 = 2, 
compared with the perturbative formulae — see text. 

In order to discuss the nonadiabatic case, where the 
polaron formation is ruled by the parameter a 2 , the re- 
sistivity data are illustrated in Fig. [3| for different val- 
ues of 7, and fixed a 2 = 10. We again recognize the 
three regimes of polaron transport (coherent, activated, 
residual scattering). The presence of such a "peak-dip" 
structure is therefore independent of the adiabaticity ra- 
tio 7 and exists whenever the carriers are of polaronic 
nature. At large 7, the resistivity obeys the following 
scaling property 

p(T,a 2 , 7 )=7/(T/ Wo ,a 2 ) (5) 

as shown in Fig. |3J Although the numerical integrals 



involved in eq. ijSJ) do not lead to an analytical expression 
for p(T), one can take advantage of the above scaling 
to derive a tractable interpolation formula valid in the 
nonadiabatic regime. To do so, we choose to separate 
arbitrarily the "coherent" part pc and the "activated" 
part pH by enforcing Matthiessen's rule p~ x = Pc +Ph ■ 
Introducing y = T/uq, the DMFT data are well described 
by 

p c (y) = A 1 a i ye a2 - 1/y (6) 
PH (y) = B7t/ 3/2 exp[A(y)/2y] (7) 

with a temperature dependent activation gap 

A(y)=a»(l-«)*^ (8) 
c/y 

and with A = 3.82, B = 4.77, c = 0.37, 5 = 0.26. 17 
Note that the form of the prefactors in eqs. ©-(0 is 
constrained by the scaling relation JSJ. The resulting 
curve for p is plotted in Fig. |3]for a 2 — 10 (full line)^ 
Let us focus on the activated regime, as was done pre- 
viously in the adiabatic case. The problem of polaron 
transport in this case is generally addressed from the 
"atomic" limit 7 — > coiiA4, treating the band parame- 
ter D as a perturbation. This yields the nonadiabatic 
textbook formulaic 

p(T) = B' 1 2 ay^ 2 e^[a 2 /2y] (9) 

valid at T > ujo, with B' = (2 7 /n) 1 / 2 [a generalization 
to lower temperatures was given in refi, eq. (97)]. Since 
the effective expansion parameter which rules the pertur- 
bative treatment is 772 itself, eq. © holds for r)2 <C 1, a 
condition opposite to eq. Q . 1419 

The inset of Fig. 3 shows an Arrhenius plot of the re- 
sistivity in the activated regime for a 2 — 10 and 7 = 2. 
As in the adiabatic marked downturn appears 

below T « wo, indicating the onset of phonon quantum 
fluctuations. We infer that this ubiquitous phenomenon 
is deeply related to the occurrence of a resistivity max- 
imum — i.e. of the very presence of small polarons — 
as it takes place both in the adiabatic and nonadiabatic 
regimes. 

In the same inset, we have also drawn the resistivity 
given by eq. © (dashed line), and its low temperature 
generalization (dotted line). Compared to the DMFT re- 
sults, we see that the perturbative formulae wildly over- 
estimate both the absolute value of the resistivity and 
the activation energy — the slope of the curve — within 
the activated regime. Besides, a closer look at eq. © 
shows that it does not obey the scaling formula (JSJ) • The 
disagreement is surprising, in view of the fact that the 
chosen parameters (772 ~ 0.01) lie well inside the range 
of validity of the perturbative approach. 

The large discrepancy comes from the narrow-band 
character of the polaronic excitation spectrum. In the 
limit D — > 0, the electron states are essentially indepen- 
dent on different sites. This, together with the fact that 
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the phonons are assumed to be local and dispersionless, 
prevents any transfer of energy between sites (the spec- 
tral function is a distribution of delta peaks), and makes 
the transition probabilities singular. Holstcin healed the 
singularity by introducing ad hoc a sizeable phonon dis- 
persion Aujph ^ 0, yielding eq. However, especially 
in narrow band materials, the optical phonons often ex- 
hibit rather weak dispersions. Obtaining a finite result 
when Au> p h — > requires to treat the electron dispersion 
(i.e. the finite bandwidth, D ^ 0) nonperturbatively. 
This can be achieved by the DMFT, as is testified by the 
finitcncss of our results^ 

In summary, we have applied the DMFT to study the 
transport properties of small polarons. The different be- 
haviors expected by standard polaron theory — coherent, 
activated and residual scattering regime — are recovered 
within a unified treatment, although notable deviations 
from the commonly accepted formulae are found. First 
of all, a broad intermediate temperature regime emerges, 
regardless of the adiabaticity parameter 7, where the re- 
sistivity is still semiconducting-like, but is strongly in- 
fluenced by phonon quantum fluctuations. This regime, 
comprised between the resistivity maximum Tf, and the 
phonon frequency ujq, should be easily detected exper- 
imentally as a downturn in the Arrhenius plots of the 



resistivity. Secondly, in the nonadiabatic regime, the 
DMFT results obey a simple scaling property, which is 
not compatible with the standard polaronic formulae of 
Holstein. Accordingly, large quantitative discrepancies 
arise in the predicted resistivity, which could result in 
wrong estimates when extracting microscopic parameters 
from the experiments. 

Finally, we would like to point out that the adiabatic 
procedure which was successfully applied to study po- 
laronic systems at finite density (half filling) within the 
DMFT scheme^ is not suitable at low densities, where 
it is not able to reproduce the correct activated behavior 
(PJ) . The reason is that the cooperative mechanism which 
leads to polaron localization at half filling — the presence 
of a bimodally distributed random field — ceases to be 
effective at low electron densities, where the phonons are 
not renormalized. In the latter case, polaron localization 
is a single-particle effect, the carriers being self-trapped 
in their own local deformation. The relative inportancc 
of the two different mechanisms of polaron trapping can 
be modulated by the density. The existence of an exper- 
imentally observable density crossover in polaronic sys- 
tems will be the subject of future investigation. 

We acknowledge M. Capone for a critical reading of 
the manuscript. 
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systems, and use the prescription given in [A. Chattopad- 
hyay et al, Phys. Rev. B 61, 10738 (2000)]. However, since 
the Kubo formula involves averages over the whole density 
of states, the results should not depend dramatically on 
the shape of N e . 

The location of the resistivity maximum (for any a 2 in the 
nonadiabatic regime) is in good agreement with the value 
given in ref.- Tt — ujq/{2 log a 2 ), while it fundamentally 
contradicts Holstein's estimate of refi . The latter is based 
on the assumption that the polaron bandwidth is reduced 
with increasing temperature, which is a drawback of the 
perturbative approach (any narrow feature in the excita- 
tion spectrum, such as the polaron band, should rather get 
broadened by thermal effects). 

Remarkably, at 0? = 10, the scaling property © holds for 
7 > 0.5 which, taking T ~ u>o, corresponds to 772 J; 1- Note, 
however, that the low temperature coherent regime is not 
ruled by eq. 0, which relies on the analysis of incoherent 
hopping processes. 

A simple calculation shows that, if A(y) is a collection of 
peaks of width Au;, the resistivity is reduced by a factor of 
Aoj/ojo relative to Holstein's result. The observed scaling 
behavior JSJ is recovered by noting that, if the broaden- 
ing of the peaks is of electronic origin, then Aw E i oc D. In 
practical cases, where both the electrons and the phonons 
disperse, it is likely that the transport properties are de- 
termined by the larger of Auj e i, Au> p h- 



